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Using a newly developed first-principles Wannier- states approach that takes into account large on-site 
Coulomb repulsion, we derive the effective low-energy interacting Hamiltonians for several prototypical high-T c 
superconducting cuprates. The material dependence is found to originate primarily from the different energy of 
the apical atom p~ state. Specifically, the general properties of the low-energy hole state, namely the Zhang-Rice 
singlet, are significantly modified by a triplet state associated with this p z state, via additional intra-sublattice 
hoppings, nearest-neighbor "super-repulsion", and other microscopic many-body processes. Possible impli- 
cations on modulation of T c , local superconducting gaps, charge distribution, hole mobility, electron-phonon 
interaction, and multi-layer effects are discussed. 
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Introduction. — The origin of high-T c superconductivity 
(HTSC) remains under fierce debate notwithstanding mon- 
umental effort for two decades (H. This fascinating phe- 
nomenon is achieved when a number of layered copper ox- 
ides are doped away from their Mott insulating parent phase. 
Although it is generally agreed that the most relevant elec- 
tron behavior is confined within the common (Cu(D 2 ) 2 ~ plane, 
ymax (j,^ ^ optimal doping) strikingly varies from 28 K 
in Ca^NazCuOzClz to 135 K in HgBa 2 Ca 2 Cu 3 8+(5 by 
modulation of the layering pattern along the less essential 
third direction. Hence, clarifying the material dependence 
of the in-plane electron behavior, especially in unbiased first- 
principles-based approaches, is an essential and effective step 
toward the resolution of the HTSC mechanism and the quest 
for higher T c superconductors. 

A recent influential advancement along the line was made 
by Andersen and coworkers yj]. Within the local-density 
approximation (LDA) of density-functional theory, they de- 
rived a one-band noninteracting Hamiltonian for several hole- 
doped cuprates, showing that —t'/t (t and t' are the first- and 
second-nearest-neighbor hopping integrals, respectively) was 
strongly material-dependent and correlated with r™ ax . Taken 
as the kinetic part of the effective one-band t-J or Hubbard 
model, the most studied model for the CuC>2 plane J^, this 
single band has been widely used to compare with angular- 
resolved photoemission spectroscopy (ARPES) 0.1110]. 

A particular puzzle thus caused is that the LDA results of 
-t'/t = 0.18 for La 2 Cu0 4 and 0.12 for Ca 2 Cu0 2 Cl 2 H 
are in striking contradiction to their diamond- and square-like 
Fermi surfaces, respectively, as observed in ARPES HHH]. In 
addition, the low-energy carriers in hole-doped cuprates are 
well accepted to be closely related to the Zhang-Rice singlet 
(ZRS) that naturally appears in deriving the t-J model from 
an interacting many -body picture |7||. This two-body singlet 
state is intrinsically different from any single-particle state. 
Therefore, it is essential to reexamine the material dependence 
by moving one step forward to derive an effective interacting 
Hamiltonian in a first-principles approach that takes into ac- 
count, from the beginning, the large Coulomb repulsion on the 



Cu sites. New material-dependent interaction effects are then 
expected to be revealed from the strong electronic interaction 
and the charge-transfer nature (H. 

In this Letter, using a recently developed first-principles 
Wannier-states (WSs) approach 

urn 

which takes into ac- 
count large on-site Coulomb repulsion, we derive the effective 
low-energy interacting Hamiltonians for several prototypical 
high-T c superconducting cuprates. The only significant mi- 
croscopic difference among the cuprates is found to originate 
from the apical atom p z state. Strong impact of the apical 
atom on t' is observed; in particular, —t'/t for Ca 2 Cu0 2 Cl 2 
is found to be considerably larger than for La 2 CuC<4, in good 
agreement with ARPES. Moreover, we show that besides the 
t' and electrostatic effects, the realistic variation of the en- 
ergy level of the apical p z states is able to tune substantially 
the local site potential and intersite "super-repulsion," directly 
modulating the local superconducting gap and charge distri- 
bution. This provides a natural explanation of recent spectro- 
scopic ima ging scanning tunneling microscopy (SI-STM) ob- 
servations ll ll Il2ll . Finally, implication on a new realization 
of the electron-phonon coupling is briefly addressed. These 
important findings should shed new light on the general mate- 
rial dependence and microscopic understanding of HTSC. 

Methods. — A three-step approach is employed to system- 
atically reduce the energy scale of the relevant Hilbert space: 
(i) The full-energy electronic structures are obtained within 
the LDA+U approximation 11311 known as a state-of-the-art 
generalization of LDA to include strong local interaction, (ii) 
At the intermediate-energy scale (~ 10 eV covering the Cu 
3<i, O 2p, and relevant apical orbitals), an effective five-band 
interacting Hamiltonian, H 5b , is derived Jgt]. The local part of 
H 5b , which includes the leading terms (on-Cu-site interaction 
and local d-p hybridization), is then diagonalized for all the 
doping levels 1 14, 15, (iii) At the low-energy scale (~ 1 
eV), an effective one-band Hamiltonian, H lb , is derived from 
canonical transformation to project out high-energy states up 
to the second order iflrjll . 

Since the apical (out-of-plane) coordination of the planar 
Cu cations is the main structural variation relevant to the 
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physics of the C11O2 plane, we present here a study of a num- 
ber of cuprates representing different apical limits: L^CuOj 
with two apical oxygen atoms per copper site, Ca 2 Cu02Cl2 
with the apical p z energy level, e Pz , moved away from the 
Fermi level by substitution of CI for O, Sr2CuC<2F2 with 
e Pz moved farther away, and Nd2Cu04 without apical atoms 
(equivalently, e Pz = 00). 

The most relevant WSs are shown in Fig. [Ha), including the 
Cu 3d, O 2p, and apical p orbitals, spanning an energy window 
of ~ 10 eV [c.f. Fig. |TJb)]. Since these WSs almost com- 
pletely decouple from other WSs, a five-band Hamiltonian 
H 5h with on-site Coulomb and exchange interactions among 
the Cu d x 2 _ y 2 and d 3z 2 _ r 2 orbitals can be unambiguously de- 
rived by matching its self-consistent Hartree-Fock expression 
with the Wannier representation of the LDA+J7 Hamiltonian 
Jgt]. The details will be presented elsewhere. Indeed, among 
the cuprates, the most noticeable difference in the parameters 



ofH 5h is e Pz 11711 . as presented in Table U 

The likelihood of a one-band picture is readily revealed by 
our WSs analysis of the LDA+U band. Fig.QJb) shows such 
an analysis for Ca2CuC<2Cl2, for there exist excellent ARPES 
data to compare with. Note that the lowest unoccupied band 
is mainly of Cu 3d x 2_ y 2 character and the highest occupied 
band is mainly of O P s character. Apparently, these two in- 
plane bands are the most relevant to the low-energy physics; 
indeed, their overall shapes agree well with the ARPES fl. 
In addition, their corresponding WSs have opposite spin char- 
acters, as clearly shown in the right panel of Fig.[TJb). Thus, 
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FIG. 1: (Color online) (a) Wannier functions centered at Cu sites 
for H 5b : O P s and P a are constructed by maximizing the weight 
of the 61- and ai -symmetrized combination of four neighboring pla- 
nar O p a orbitals, respectively; P z is the a\ -symmetrized combina- 
tion of two apical p~ orbitals. (b) Left panel: Band structures of 
Ca2Cu02Cl2 from LDA+U calculations (dots) and Wannier states 
analysis (lines). r=(0, 0, 0), X=(w, 0, 0), M=(f , § , 0). Right 
panel: Corresponding local energy level splitting due to interactions 
and d-p hybridization, filled with electrons (solid arrows) or holes 
(dashed arrows). Pz are derived from two highest a\ hybrids. 



the simplest one-band picture is to view these two anti-spin- 
polarized bands (both are antibonding between Cu 3d x 2_ y 2 
and O P s WSs) as the effective upper and lower Hubbard 
bands (EUHB and ELHB) with the charge-transfer gap as the 
effective Hubbard gap. 

A more appropriate approach is to project out high-energy 
states systematically, based on exact diagonalization of the 
leading terms (here the local part of H 5b ) 1 14 , 15 , 13] . Specif- 
ically, our approach benefits greatly from having all the^zraf- 
principles WSs, including the O 2p states, centered at the Cu 
sites [see Fig. 02 a )] by construction. This reduces the deco- 
rated CuC<2 lattice to a simple square lattice, and unambigu- 
ously defines the local Hamiltonian to be diagonalized. 

It is appropriate to keep in the low-energy space only the 
local one-hole and two-hole ground states, referred to as \a) 
and |ZRS), respectively, as far as the mechanism of HTSC is 
concerned, for which the relevant energy scale is expected to 
be smaller than 0.5 eV 13]. Note that the first-excited two-hole 
state is a spin-triplet (referred to as |axial)) composed of one 
61 hole and one a\ hole with excitation energy A2 = 0.81 eV 
for Ca2CuC<2Cl2, for example. A one-band model is thus ob- 
tained by projecting out all the other states by canonical per- 
turbation and then mapping to a constrained fermion system 
without double occupancy llal . 

Additional caution was paid to La2CuC<4. With a model 



tetragonal structure 



m was paic 

QElIElI, 



we obtain e 



Pz 



'Ps 



-0.52 



eV and A2 is merely 0.13 eV, suggesting that a two-band 
Hamiltonian including |axial) be more appropriate. For this 
subtle case, other more realistic considerations are important; 
for example, we find that the 5° octahedral-tilting in the real 
material increases e Pz — e Ps by 0.25 eV and A2 by 100%. 
Such a tilting effect is interesting, as the oxygen dopants in 
Bi2Sr2CaCu208+,5 also induce the tilting of nearby oxygen 
octahedra. To account for such effects, we present in Table U 
the results for the model La2CuC<4 with e Pz — e Ps being sim- 
ply increased by 0.8 eV (now A2 = 0.5 eV); the e Pz effects 
are still large though considerably abated. 

Results and Discussion. — The resulting effective one-band 
t- J-like Hamiltonian is given by lEoll 



H 



lb 



J2^(cl^ + H.c.) + jJ2(S l -S J 



) 



1 



(ijk)tr 

(ij) i ij 



(1) 



where c ity is the annihilation operator for the constrained 
fermion with spin er at site i, — <? ia c ir! , rii = n ia, 

and Si = ch&nvCiff with °W bein § the Pauli matri- 
ces. All the site pairs (ij) stand for nearest neighbors ex- 
cept tij extends to the third nearest neighbors (referred to as 
t"). Note that the first two lines of Eq. (Q]i resemble the t- 
J model mapped out from the one-band Hubbard model and 
have been extensively studied 12 ill . The other terms include 
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s Pz -dependent "super-repulsion" Vij, site potential p^ and 
energy "constant" . 

The derived parameters for several prototypical cuprates 
upon hole doping are listed in Table U For comparison, the 
contributions purely from the b\ orbitals (from which ZRS 
is formed) are given inside the parentheses in Table Q] Re- 
markably, there exists practically no material dependence if 
only the b\ orbitals are considered. The universality for t and 
J survives inclusion of the a\ orbitals, since they are solely 
determined by the b\ orbitals, suggesting that the magnetic 
mechanism be responsible for the robustness of superconduc- 
tivity among the cuprates. 

More amazingly, the material dependence of these materi- 
als almost entirely results from the influence of the apical P z 
orbitals; the other a\ orbitals have negligible contributions 
as made clear from comparing all the cases with Nd2Cu04. 
t' , t", Vij, pi, and £ij all show strong material- (s Pz -) de- 
pendence [c.f. Fig. [21a)]. First of all, three-site kinematical 
processes via virtual nearest-neighbor hoppings to the P z or- 
bital dramatically renormalize t' and t", as suggested in Ref. 
1 1411 . These modulations would seriously affect the mobil- 
ity of hole quasiparticles and the shape of the Fermi surface, 
given that the nearest-neighbor hoppings of hole quasiparti- 
cles are suppressed by the antiferromagnetic spin correlation 
1 2111 . Note that —t'/tis substantially smaller in La2Cu04 than 
in Ca2CuC>2Cl2, consistent with their observed diamond- and 
square-like Fermi surfaces, respectively flllt]. 

Quantitatively, our calculated values of — t'/t for the above 
two compounds are considerably smaller than the widely used 
ones (~ 0.15 and 0.3, respectively) estimated from fitting 
ARPES in the very lightly doped t-J model fl H HQ. 
Interestingly, the single-hole quasiparticle dispersions calcu- 
lated with the present parameters for A2CUO2X2 {A=Ca, Sr; 
X=F, CI) in the self-consistent Born approximation (SCBA) 
10 still agree well with available ARPES data fH, as shown 
in Fig. Ob). This is because the dispersion is actually more 
sensitive to t"/t, which is consistent between ours and pre- 
vious estimations. Moreover, the fitted values of t' and t" 



were shown to depend sensitively on inclusion of the three- 
site hopping terms oc J/4 J22I1 . Our first-principles derived 
Hamiltonian finally provides an unambiguous benchmark of 
these material-dependent parameters. 

In sharp contrast, the "downfolded" LDA one-band led to 
a puzzling opposite trend {—t'/t = 0.18 for La2CuC<4 and 
0.12 for Ca2CuC<2Cl2) \M that appears in contradiction with 
ARPES and our parameters. This is, however, not alarming, 
since i's in Eq. (1) describes hoppings of ZRS (intersite swap- 
ping of the |ZRS) and \a) states Jj), while the LDA t's de- 
scribes hopping of the single-particle LDA \a) state. This in- 
trinsic difference can be easily observed from Fig.[TJb) where 
EUHB and ELHB are shown to have different WSs characters, 
an essential feature missed in the picture of the LDA one- 
band plus effective Hubbard interaction ]2[]. Clearly, build- 
ing the strong many -body characteristics in first-principles ap- 
proaches before the downfolding HHQjl], as presented here, 
is necessary for the low-energy physics of cuprates. 

Next, the striking appearance of the unconventional "super- 
repulsion" Vij originates from two-site kinematical processes 
as demonstrated in Fig. [3] For a purely one-band system in the 
large on-site repulsion limit [Fig. a)], the virtual hoppings 
between two singly-occupied neighboring sites generate the 
well-known superexchange effect —jn ia ?ij S . Given an extra 
fully-occupied orbital to be projected out (P z in the present 
work), additional virtual kinematical processes give rise to 
three new spin-independent terms [oc voq, v±o, and v\\, see 
Fig. [3b)-(d)] to the targeted one-band Hamiltonian, leading 
to Vij = 2v w - voo - uii, pt = 4(v o - v w ) and e ?J = -v o- 

This new effective repulsion, V^ (named "super-repulsion" 
in analog to superexchange), is to be distinguished from di- 
rect Coulomb interaction: Not only is Vij controlled by e Pz , 
but also its virtual kinematical origin makes it less subject to 
the electronic screening. Besides, as presented in Fig. EJc), 
Vij rapidly increases as \t' /t\ decreases from Nd2CuC»2 to 
La2Cu02. Therefore, for any realistic study of the cuprates, 
the Vij effect should be addressed along with the t' effect. 

As for superconductivity, the repulsive Vij apparently 



TABLE I: Derived parameters of Eq. (TJ for four prototypical 
cuprates upon hole doping. Inside the parentheses are contributions 



purely from the 61 orbitals 


The energy 


unit is meV. 






Nd 2 Cu0 4 


Sr 2 Cu0 2 F 2 


Ca 2 Cu0 2 Cl 2 


La 2 Cu0 4 


ymax (K) 


N/A 


46 


28 


38 


e Pz ~ e Ps 


00 


3010 


810 


280" 


Jlda+c/' 


131 


155 


145 


144 


J/2 


138 (138) 


138 (138) 


131 (131) 


129 (129) 


t 


431 (431) 


467 (467) 


459 (459) 


488 (488) 


t'J\t\ 


-0.33 (-0.35) 


-0.25 (-0.32) 


-0.19 (-0.32) 


0.01 (-0.30) 


t"/\t\ 


0.23 (0.24) 


0.19(0.23) 


0.16 (0.23) 


0.06 (0.22) 


Vij 


29 (27) 


39 (29) 


54 (30) 


137 (31) 




-795 (-800) 


-337 (-394) 


-295 (-412) 


-354 (-529) 


Sij 


-36 (-14) 


-86 (-17) 


-133 (-17) 


-258 (-19) 



"Increased by 0.8 eV (see text). 

*Estimated from total energies of the (anti)ferromagnetic states. 
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FIG. 2: (a) Parameter changes by the cii orbitals as a function of 
e Pz . (b) Calculated one-hole quasiparticle dispersions compared 
with ARPES on Sr 2 Cu0 2 Cl 2 (open circles) H. (c) Vij versus \t' /t\. 
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weakens the local pairing strength, regardless of the actual 
pairing mechanism, given that paired holes tend to reside as 
nearest neighbors (implied by the d x 2_ y 2-wave symmetry of 
the order parameter [23]). This provides a natural avenue 
for modulating the local superconducting gap ll ill as well as 
T c . On the other hand, the suggested trend of T c max ~ t'/t 
II2L |l4Jl — which turned out to be controversial in several nu- 



reveal a strong correlation between the local superconducting 
gap and Cu-O ap icai distance 112 1. 



merical studies 124 I25L 12611 — is clearly violated by compar- 
ing Caa-^Na^CuOaClz (T c max =28 K) with La^Sr^CuCU 
(T c max =38 K). It is thus fair to speculate that the full material 
dependence of T™ ax involves other complex aspects of real 
materials (e.g., sample cleanness, doping approaches, two- 
band scenarios). Further investigation is needed to identify 
those aspects that affect T" lax beyond our present discovery. 

With regard to charge distribution (CD), modulation of \ii 
tends to induce charge inhomogeneity. As shown in Fig. [2a), 
Hi depends strongly on the apical environment. This explains 
strong interlayer charge inhomogeneity in multilayer systems 
such as HgBa2Ca n _iCu„02n+2+(5 with n > 3 1127 [and atom- 
ically perfect Lai.85Sro.i5CuC»4/La2Cu04 films 12811 . as [ii in 
distinct CUO2 layers can be considerably different. On the 
other hand, repulsive Vij, which disfavors hole accumulation, 
effectively counters /i; to give a smoother in-plane CD. 

The present results would shed new light on the understand- 
ing of other exotic features of cuprates. For example, recent 
SI-STM experiments on I^S^CaC^Os+i indicated that lo- 
cal superconducting gaps vary from 20 — 70 meV correlat- 
ing with oxygen dopants, Os, while low-energy charge den- 
sity variations are weak 111 111 . These striking results have been 
phenomenologically attributed to a strong local modulation 
of electron pairing strengt h |2 9ll. J in the t-J model ll3(il . or 
Hi with two types of O5 13 ill . The present studies suggest 
that Og perturb O ap i ca i and thus significantly modulate Vij, 
instead of J, in addition to fj,i, t 1 , and t" . This new scenario 
also agrees well with more recent SI-STM measurements that 
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FIG. 3: Schematics of virtual kinematical processes and correspond- 
ing effective one-band interaction terms for (a) Heisenberg superex- 
change within a one-orbital system and (b)-(d) additional "vacuum 
fluctuations" in a two-orbital system. Electrons (solid arrows) in the 
targeted and projected out orbitals are black and gray, respectively. 
Dashed lines denote hopping paths. 



Furthermore, the strong e Pz dependence of t' , t", Vij, and 
Hi implies the electron-phonon interactions involving vibra- 
tion of the apical atoms may be strong. In particular, the 
V^ variation in this way yields an unconventional higher- 
order electron-phonon interaction. The determination of their 
actual strengths will be pursued elsewhere. The enhanced 
electron-phonon interactions could cooperate with Vij in fa- 
voring other competing orders such as "stripe," especially in 
doped La2Cu04 13211 where we have shown that the pair- 
breaking Vij is comparable to J. 

Finally, our findings point to the imp ortance of examining 
HTSC in multi-layer systems 1271 13 311 or cuprate superlat- 
tices 128ll . as the influence of the apical P z orbitals will be 
significantly layer-dependent and can be controlled (e.g., via 
strain and field effects) to tune hole concentration, mobility, 
and super-repulsion among the CUO2 layers, a very interest- 
ing problem for future investigation. 
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